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ABSTRACT 


A  flat  plate  experiences  a  transverse  acceleration  under  the  action 
of  a  pressure  that  also  produces  plastic  yielding  throughout  the  plate. 
Taylor  instability  is  considered  on  the  assumptions  of  a  single  egree  of 
freedom,  incompressible  flow,  and  the  Prandtl-Reuss  constitutive  equations 
in  the  limiting  cases  of  (A)  infinitesimal  perturbations  and  (B)  rigid- 
plastic  flow  (yield  stress/shear  modulus  -♦  0). 
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1.  Introduction 

A  flat  plate  is  subjected  to  a  uniform  pressure,  say  p  ,  over  its 
lower  surface,  say  y  =  0,  and  is  otherwise  unconstrained.  The  resulting 
acceleration, 

g  =  P0M  >  (1.1) 

is  directed  along  the  positive  y-axis  in  a  reference  frame  that  moves  with 
the  plate;  p  is  the  mean  equilibrium  density,  and  h  is  the  equilibrium 
thickness.  We  consider  the  dynamic  stability  of  this  equilibrium  configura¬ 
tion  with  respect  to  a  small ,  periodic  perturbation,  such  that  the  displace¬ 
ment  of  the  lower  face  is 

\(x,t)  =  q(t )  coskx  (|  q|  «  h)  (1.2) 

in  a  reference  frame  that  moves  with  the  plate. 

We  assume  that 

Po»s1,  (1.3) 

where  s^  is  the  yield  stress  in  shear,  and  that  the  material  obeys  the 
Prandtl-Reuss  constitutive  equations,  with  the  yield  stress  and  shear 
modulus,  G,  as  constitutive  parameters.  We  also  assume  that  the  material 
is  homogeneous  (thereby  neglecting  the  y-dependence  of  the  equilibrium 
density  p)  and  incompressible.  Admitting  only  a  single  degree  of  freedom 
and  invoking  energy  considerations,  we  establish  the  differential  equations 
for  q(t)  in  two  limiting  cases,  namely: 

(A)  (G/SjMql^-0 

and 


(B)  s1/G  -  0  . 


2.  Kinematics 


Taylor's  resulos  for  an  accelerated  slab  of  inviscid,  incompressible 
fluidx  yield  the  normal  modes 

$i(x,y,t)  =  f±(t)e±kycoskx  (2.1) 

for  the  velocity  potential  corresponding  to  a  disturbance  of  wavenumber  k. 

The  $  mode  is  unstable,  with  an  exponential  growth  rate  of  (gk)^j  the  $ 
mode  is  stable,  with  an  angular  frequency  (gk)^.  Guided  by  these  results,  w 
we  base  our  analysis  on  the  single-degree-of-freedom  approximation 

§(x,y, t )  =  -k_14(t )e"kyc.oskx  ,  (2.2) 

which  satisfies  the  boundary  condition  implied  by  (1.2),  namely 

d§/3y  =  3T|o/dt  (y  =  T]o  =  0)  .  (2.3) 

The  corresponding  motion  of  the  upper  surface,  y  =  h,  is  given  by 

yx,t)  =  .  (S.k) 

The  potential  energy  associated  with  the  displacements  of  (1.2)  and 
(2.1+)  in  the  accelerated  reference  frame  is  given  by 

dV/dx  =  fcpg(T^  .  T®)  (2.5a) 

=  -  £pg(l  "  e"2kil)q‘_(t)cos2kx.  (2.5b) 

Averaging  (2.5b)  over  one  or  more  wavelengths,  we  obtain 

<V>  =  -  ipg(l  -  e"21dl)q2(t)  .  (2.6) 

The  corresponding  average  of  the  kinetic  energy  is 

r k  o  p 

<£>  =  h  I  <  +  r  >  dy  (2,.7a) 

=  4Pk-1(l  -  e"2kh)42(t)  .  (2.7b) 

G.I.  Taylor,  Proc.  Roy.  Soc.  (London)  A201,  192(1950);  Scientific 
Papers  (Cambridge  University  Press,  i960)  vol.  3,  p.  532.  Note  that  the 
acceleration  in  Taylor's  paper  is  directed  along  the  negative  y-axis,  so  that 
the  roles  of  in  respect  to  stability  are  reversed. 
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3.  Constitutive  Equations 

Let  s^  and  e^^  be  the  stress  and  strain-rate  deviators, 

J2  =  SiJSij  and  Z2  =  ^ij^ij 
their  second  principal  invariants,  and 


(3.1a,b) 


W  = 


s, ,e. , 
i<3 


(3.2) 


the  rate  at  which  work  is  done  in  deforming  a  unit  volume  of  material.^  We 
recall  that,  by  definition,  the  first  principal  invariants  vanish; 


The  Prandtl-Reuss  hypotheses  imply  J  £  2s 
yield  condition.  The  constitutive  equations  are^ 


=  0  . 

2 

^  ,  with  J 


2 


(3.3a,b) 


as  the 


s, ,  +  =  2Ge.  . 

ij  ij 


where 


(3.4) 


X  =  (G/s2)W  (J2  =  2s2  &  W>  0)  (3.5a) 

=0  (J2  <  2s2  or  J2  -  2~2  &  W  <  0)  .  (3.5b) 

2  •  2 
The  conditions  J2  =  2s^  and  W  >  0  imply  plastic  flow;  J2  <  2s^  implies 

elastic  deformation;  J2  =  2s2  and  W  <  0  imply  unloading  at  the  yield  limit. 

The  Prandtl-Reuss  equations  are  indeterminate  for  a  state  of  static 
equilibrium  on  the  yield  surface.  In  the  present  instance,  the  plastic 
compression  of  a  flat  plate,  we  invoke  the  yield  condition  and  symmetry 
considerations  to  obtain 


s  =  s  =  sn  4/3  ,  s  =  -2s  .4/3  ,  s  =  s  =  s  =0  (3.6a) 

xx  zz  l/w  ’  yy  ’  xy  yz  xz  '  ' 


We  note  that  the  development  in  §3,  and  also  in  $4  through  (4.9b), 
does  not  require  the  assumption  of  incompressibility.  The  invocation  of  this 
assumption  renders  the  total  strain-rate  tensor  and  the  strain-rate  deviator 
identical. 

2 

W.  Prager  &  P.  G.  Hodge,  Theory  of  Perfectly  Plastic  Solids  (John 
Wiley  8s  Sons,  New  York,  1951 ),  pp.  27-30. 
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or,  more  briefly, 

sij  =  (si/^>6u  }>  -2'  ^  >  (3-6b) 

where  1,  -2,  1^  comprises  the  non-zero  terms  in  a  diagonal  tensor.  We 
note  that  these  results  overlook  a  small  neighborhood  of  the  upper  surface, 
of  the  order  of  (s1/pQ)h  in  thickness,  in  which  the  material  has  not  yielded. 
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4.  Small -perturbation  Approximation 

We  now  consider  small  perturbations  about  the  equilibrium  state  of 
(3.6),  say 


•  s<°>  +  .W  (l  s<1>1  «  .  1 

-i)  -  Sij  +  Sij  |Sij  1  ^  Sl|  ’ 


(4.1) 


where  s)^  is  the  perturbation  stress  associated  with  e44.  The  first 

1J  #  lj  "" 

approximation  to  W  on  this  assumption  is 

ft'1*  -  s^h 

w  -  su  eij 

=  (s/V3)(e  -  2e  +  e  ) 

'  r  ,v  xx  vy  zz' 


(4.2a) 


(4.2b) 


-  -  </3s1e 


(4.2c) 


where  (4.2b)  and  (4.2c)  follow  from  (4.2a)  by  virtue  of  (3.6b)  and  (3.3b). 
The  corresponding  first  (linearized)  approximation  to  (3.4)  is 


sJV  +  -  2Ge  , 


(4.3) 


where  the  square  bracket  imply  the  positive  part  of  the  bracketed  term--i.e. 


[W]  =  W  (W>  0) 

=  0  (W  <  0) 

Substituting  (3.6b)  and  (4.2c)  into  (4.3),  we  obtain 

sf^  =  2Ge.  -  G[-e  ]6.,  SI,  -2,  1,  . 
ij  ij  yy  ij  ^  ’  J 

Integrating  (4.5,  we  obtain 


(4.4a) 

(4.4b) 


(4.5) 


■S>  -  2Gei.)  -  o6u 


t1’  -2-  b  J  [-4y]4t  * 

O 


(4.6) 


where  we  define  e . .  and  t  such  that 

ij  o 


s  ;  i  =  e *  ,  =  0  at  t  =  t 
lj  IJ  o 


(4.7) 
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We  require,  for  the  construction  of  the  linearized  equations  of 
motion,  the  second  (quadratic)  approximation  to  W,  say 

W  =  +  w^2)  . 

Multiplying  (4.6)  through  by  e^  and  invoking  (3.3b),  we  obtain 


(*.8) 


W^  =  2Ge,  ,e,,  -  G(e  -  2e  +e  )  f  [_e  ]dt 
ij  ij  xx  yy  zzM  yyJ 

o 


(4.9a) 


=  2Ge  e  ,  +  3Ge  /  [-e  ]dt  . 

ij  yy  yy 


(4.9b) 


The  strain-rate  tensor  implied  by  (2.2)  is 

-  4,1J  e  k4(t)e-kyMiJ(x)  , 


where 


(4.10) 


coskx  -sinkx  0 
=  -sinkx  -coskx  0 
0  0  0 

Substituting  (4.10),  together  with  its  temporal  integral  from  tQ,  into 
(4.6)  and  (4.9b),  we  obtain 


1  ^  =  0ke_ky  2q(t)H1j  -  6i J  |1,  -2,  lj-  f  [4('t)ooskx]dt  (4.12) 

'  I 


W^  =  Gk2e"2kyq(t) | 4q(t )  -  3coskx  J  [q(t)  coskx]dt  ,  (4.13) 


where,  by  definition, 


q(t0)  =  0  . 


(4.14) 
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Comparing  (4.12)  and  (4.1),  we  Infer  that  (4.12)  and  (4.13)  are  valid 
approximations  if 


k|  q|  «  s  /g 
1  1  max  ]/  u  > 


(^.15) 


as  anticipated  in  1. 


Remarking  that  the  mean  value  of  vani  qVioq  -k  j 

.  01  w  vanishes  by  virtue  of  the 


periodicity  of  e  ,  we  obtain 


<ii>  =  <w(2)  >  =  |  Ok2e-%(t)4(t) 


('1.16) 


Substituting  (4.16),  together  with  (2.6)  and  (2  7b)  in+o  tb* 

equation  V  '  ^.7 bj,  into  the  energy-balance 


a  .  rh  . 

3t  ^  +  )  +  I  <W>dy  =  0 

^  o 

and  dividing  the  result  through  by  4(p/k)(i  .  e'2kh)4,  Ve  obtain  the 
linearized  equation  of  motion 


(4.17) 


3(t)  +  ip-^k2  -  gk]  q(t)  =  0  . 


(4.18) 


Substituting  g  from  (1.1)  into  (4.18),  we  infer  the  stability  of  small 
perturbations  for  which 


44  >  6po/l3G  . 


Combining  (4.15)  and  (4.19),  we  obtain 

lqlmax/h  « 


(4.19) 


(4.20) 


as  a  necessary  condition  for  the  validity  of  the  stability  criterion  (4.10) 

This  is  a  rather  severe  rPe-i- _  i  ^  ' 

vere  restriction  in  conseauenop  n-p  +->.0  D  _ j  _ 


(1.3). 


consequence  of  the  a  priori  condition 


10 


5.  Rigid-plastic  Approximation 

We  can  avoid  the  restriction  (4.20)  by  letting  s,/G  -*  0  at  the  out- 
set  to  obtain  the  rigid -plastic  (von  Mises)  constitutive  equation?"^ 


Sij  “  s1^2/I2^  ej 


(5.1) 


W-Si(2I2)*  (5.2) 

in  place  of  (3.4)  and  (3.5).  We  emphasize  that  (5.1)  does  not  allow  for 
elastic  unloading  from  the  yield  surface  and  therefore  may  imply  stress 
discontinuities  in  both  space  and  time.  We  also  emphasize  that  W  ,  as 
given  by  (5.2),  is  non-negative. 

Evaluating  I ^  from  (3.1b)  and  (4.10),  we  obtain 

W  =  2S;Lke-kylq(t)|  .  (5.3) 

Substituting  (5.3),  together  with  (2.6)  and  (2.7b),  into  (4.17),  we  obtain 


cj(t )  -  gkq(t)  +  4(1  +  e"kh)"1(ks1/p)sgnq(t)  =  0 


(5.4) 


The  differential  equation  (5.4)  is  similar  to  that  for  a  Coulomb- 
damped  oscillator  (but  with  the  important  difference  that  the  square  of  the 
angular  frequency,  -gk,  is  now  negative).  Introducing  the  dimensionless 
variables  f  and  t  according  to 


q(t)  =  6hf (t ) 


(5.5) 


t  =  (gk)^t 

where  [s^/pgh  s  s^p  by  virtue  of  (1.1)] 


6  =  4(1  +  , 


(5.6) 


(5.7) 


'’Prager,  op.  cit.  ante,  p.  31*. 


****'■' 


we  rewrite  (5 A)  in  the  normalized  form  (dots  now  imply  differentiation 
with  respect  to  t) 
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f (t)  f (t )  +  sgnf  =  0  .  (5.8) 

We  emphasize  that  the  wave  number  no  longer  appears  explicitly  in  (5.8)  and 
that  it  enters  primarily  through  the  time  scale  (gk)"^j  the  variation  of  the 
factor  (l  +  e-^)”^,  from  ^  to  1  as  kh  varies  from  0  to  00 ,  is  small  compared 
with  the  uncertainties  implied  by  the  antecedent  approximations. 

Transforming  (5*8)  to 


df  f  -  sgnf 
df  “  -p 


we  obtain  the  integral  curves 

(f  -  sgnf)2  -  f2  =  (fQ  -  sgnfQ)2  -  f  ,  (5.10) 

where  f  and  f  are  initial  values.  These  are  sections  of  rectangular 
o  o 

hyperbolae  as  sketched  in  Fig.  2;  the  arrows  designate  the  direction  in 
which  a  given  integral  curve  must  be  transversed  (with  increasing  t). 

There  are  six  distinct  possibilities,  namely: 

(A)  fo  >  0 

(All  f  <  1  ..  f 
v  '  o  o 

(A2)  1  -  f  <f  <  o  +  f 

0  0  o 

(A3)  f0>l  +  i0i 

• 

(B)  f  <0  and 

o 

(Bl)  f  >  -1  -  f 

(B2)  -1  +  f0  <  f0  <  -1  -  f0 
(B3)  fc  <  -1  +  f„  . 


f 
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Trajectories  that  originate  in  either  (A^)  or  (B^),  i.e.  in  the  corridor 

(f  +  f)sgnf  *  1  ,  (5-11) 

must  terminate  on  the  interval  |f|  <1  and  f  =  0  and  correspond  to  stable 
motions,  whereas  trajectories  originating  outside  of  this  corridor  diverge 
according  to 

•  • 

f (t )  =  sgnfQ(l  -  coshT)  +  fQcoshT  +  f^sinhT  (5.12a) 


^(f  +  f  -  sgnf  )eT  . 

*'  o  o  o' 


(5.12b) 


Considering  further  only  those  disturbances  that  originate  from  a 
state  of  rest  (f  =  Ob),  we  conclude  that  disturbances  for  which 


q  I  >  6h  =  (q  ) 
o'  'o'er 


(5.13) 


will  diverge  according  to 

q(t)~£[|qj  -  (qQ)cr]expL(gk)^t].  (5.1*0 

We  observe  that  disturbances  of  short  wavelength  (large  k)  have  a  higher 
threshold  of  instability  (although  this  is  a  relatively  small  effect)  but  grow 
more  rapidly  if  this  threshold  is  exceeded.  We  emphasize,  however,  that  the 
basic  hypotheses  on  which  these  conclusions  are  based  cannot  remain  valid 
for  disturbances  of  very  short  wavelength;  in  particular,  viscous  effects 
must  be  stabilizing  for  sufficiently  large  k. 
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6.  Two-degree-of -freedom  Problem 

We  consider  briefly  the  extension  of  the  preceding  analysis  to 
accomodate  two  degrees  of  freedom,  with  (1.2)  replaced  by 

T[o(x,t)  =  q1(t)coskx  +  ^(tjsinkx  .  (6.1) 

We  obtain  the  corresponding  generalizations  of  <\£>,  <T>,  and  <WV  '>  simply  by 
replacing  q  by  q1  +  ^  in  (2.6),  q  by  q±  +  in  (2.7b),  and  qq  by  q^^  q^ 
in  (4.16),  in  consequence  of  which  each  of  q^  and  q^  satisfies  (4.18)  indepen¬ 
dently— i.e.,  q^  and  q^  are  normal  coordinates  in  the  small-perturbation 
approximation . 

The  position  is  otherwise  for  the  rigid-, 'las tic  approximation,  in 
that  q^  and  q^  are  coupled  through  the  plastic  forces.  The  above  generaliza- 
tions  of  <V>  and  <I>  remain  valid,  but  |q|  must  be  replaced  by  the  non¬ 
negative  radical  (q^  +  q^)^  in  (5*3)  to  generalize  W.  The  analogy  with  the 
Coulomb-damped  oscillator  continues  to  hold,  however,  by  virtue  of  the  fact 
that  the  rigid-plastic  constitutive  equation  (5*1)  postulates  a  stress 
tensor  s  that  is  parallel  to  the  strain-rate  tensor  e  but  independent  of 
its  magnitude  l|*  this  is  a  simple  generalization  of  the  constitutive 
equation  for  Coulomb  friction,  which  postulates  a  force  vector  that  is 
parallel  to  the  velocity  vector  but  independent  of  its  magnitude.  Invoking 
this  analogy,  we  find  that  the  two-degree-of -freedom  generalization  of  (5-4) 
is 

-  gkqi  +  4(1  +  e"kh)‘1(ks1/p)(q^  +  q^)-^  q±  =  0  (i  =  1,2)  .  (6.2) 

Introducing  the  dimensionless  variables  f(T)  and  0(t)  according  to 

q^  =  6hfcos9  ,  q^  =  6hfsin9  ,  (6.3) 

where  t  and  5  are  defined  by  (5*6)  and  (5.7)>  we  obtain 


••  •  p  •  p  p  #  p  — A- • 

f  -  f  -  f0  +  (f  +  f  0  )  *f  =  0 

(6.4a) 

f9  +  2fS  +  (f  +  f  0  )  7f0  =  0  . 

(6.4b) 

r 
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Introducing 


u(f)  =  f,  v(f )  =  f20. 

(6.5) 

we  obtain  the  first-order  differential  equations 

du/df  =  fu-1  +  f"3u_1v2  -  f(f2u2  +  v2)"^ 

(6.6a) 

and 

dv/df  =  -fvu-1(f2u2  +  v2)"^. 

(6.6b) 

We  may  assume  u  2:  0  without  loss  of  generality  (cf.  Fig.  2,  which 
illustrates  the  symmetry  between  motions  with  positive  and  negative  values  of 
fQ).  We  therefore  may  infer  from  (6.6a)  that  du/df  is  positive-definite,  and 
hence  that  integral  curves  in  a  u,f -plane  diverge,  for  f  >  1;  accordingly, 

<q10  +  >0ll> 

(6.7) 

a  simple  generalization  of  (5.I3),  is  a  sufficient  condition  for  instability. 
We  have  been  unable  to  obtain  a  generalization  of  (5. 11),  but  the  form  of 
the  differential  equations  strongly  suggests  that  motions  for  which  |v|  >  0 
(two  degrees  of  freedom)  are  less  stable  than  comparable  motions  for  which 
v  =  0  (one  degree  of  freedom).  On  the  other  hand,  we  infer  from  (6.4b)  that 
v  must  be  a  monotonically  decreasing  function  of  t  (unless  v  =  0),  in 
consequence  of  which  the  spatial  phase  difference  between  the  two  degrees 
of  freedom  can  have  only  a  transient  effect. 
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Fig.  2--Typical  integral  curves  of  (5.9)*  Trajectories  that  originate  within  the 
corridor  bounded  by  the  cross-hatched  lines  are  stable;  trajectories  that  originate 
outside  of  this  corridor  are  unstable. 


